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ABSTRACT 

We present the results of a systematic study using ASCA of the ionization 
state for six "mixed-morphology" supernova remnants (MM SNRs): IC 443, 
W49B, W28, W44, 3C391, and Kes 27. MM SNRs show centrally filled thermal 
X-ray emission, which contrasts to shell-like radio morphology, a set of character- 
istics at odds with the standard model of SNR evolution (e.g., the Sedov model). 
We have therefore studied the evolution of the MM SNRs from the ionization 
conditions inferred from the X-ray spectra, independent of X-ray morphology. 
We find highly ionized plasmas approaching ionization equilibrium in all the 
MM SNRs. The degree of ionization is systematically higher than the plasma 
usually seen in shell-like SNRs. Radial temperature gradients are also observed in 
five remnants, with cooler plasma toward the limb. In IC 443 and W49B, we find 
a plasma structure consistent with shell-like SNRs, suggesting that at least some 
MM SNRs have experienced similar evolution to shell-like SNRs. In addition to 
the results above, we have discovered an "overionized" ionization state in W49B, 
in addition to that previously found in IC 443. Thermal conduction can cause the 
hot interior plasma to become overionized by reducing the temperature and den- 
sity gradients, leading to an interior density increase and temperature decrease. 
Therefore, we suggest that the "center-filled" X-ray morphology develops as the 
result of thermal conduction, and should arise in all SNRs. This is consistent 
with the results that MM SNRs are near collisional ionization equilibrium since 
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the conduction timescale is roughly similar to the ionization timescale. Hence, 
we conclude that MM SNRs are those that have evolved over ~10 4 yr. We call 
this phase as the "conduction phase." 

Subject headings: conduction — plasmas — radiation mechanisms: thermal - 
supernova remnants — X-rays: ISM 

1. Introduction 

The most prominent characteristics of the plasma in a supernova remnant (SNR) are 
its high electron temperature (~10 7 K) and low density (<1 cm -3 ). It takes a substantial 
fraction of the life of a SNR (10 4 ~ 5 yr) for ions in such an optically-thin thermal plasma 
to reach collisional ionization equilibrium (CIE) with the electrons heated by the shock 
wave. Such an "underionized" state is often referred to as non-equilibrium ionization (NEI) 
and has been observed in many SNRs. Thus, investigating the range of ionization states 
in a SNR provides significant information about its evolution. The objects of the present 
study are "mixed-morphology" SNRs (MM SNRs) whose formation mechanism is not yet 
well understood. The properties of MM SNRs have been most succinctly described in Rho 
& Petre (1998). They have a shell-like radio morphology, and their steep radio spectrum 
is characteristic of synchrotron radiation from electrons accelerated by the shock wave. In 
contrast, the X-ray morphology of MM SNRs is centrally peaked, with little or no limb 
brightening. Unlike the centrally peaked synchrotron radiation from pulsar wind nebulae 
such as the Crab Nebula, the X-ray emission is thermal, showing line emission from abundant 
metals such as Mg, Si, S, and Fe. In general, thermal X-rays are thought to arise mainly 
from the shock-heated and compressed material behind the forward blast wave that would 
give rise to a shell-like morphology. Thus, the properties of MM SNRs are inconsistent with 
the initial ~10 4 yr of the standard model of SNR evolution (i.e., the Sedov model), when the 
shock front is moving fast enough to heat ambient material to X-ray emitting temperatures. 
However, MM SNRs represent about 20% of the X-ray detected Galactic SNRs, suggesting 
that the conditions for creating such remnants are not rare. A large number of them are 
interacting with molecular clouds and/or H I clouds. For example, OH maser emission at 
1720 MHz is detected from seven MM SNRs, a remarkably high fraction (Yusef-Zadeh et al. 
2003). 

Two scenarios have been invoked to explain the X-ray morphology of MM SNRs. One is 
that the interior X-ray emission arises from the gas evaporated from shocked clouds (White 
& Long 1991), implying that MM SNR have evolved under special conditions. The other 
is that as a SNR evolves, the temperature and density of the hot interior plasma gradually 
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becomes uniform through thermal conduction, and once the remnant has entered the "shell 
forming" stage in which the forward shock no longer heats ambient material to X-ray emitting 
temperatures, it becomes visible as centrally brightened X-ray emission (Cox et al. 1999; 
Shelton et al. 1999). However, the formation mechanisms of MM SNRs are still controversial. 
We investigate the ionization degree of the plasma to study the evolution of MM SNR since 
the ionization state of a remnant can yield an age estimate that is independent of the X-ray 
morphology. 

In this paper, we first describe our selection of targets for analysis of ionization state 
and plasma structure (section 2). Section 3 presents the analysis and results of individual 
remnants, the details of which are in Kawasaki (2003). In section 4, we summarize the 
results of the MM SNRs through comparison with shell-like SNRs. We then discuss possible 
mechanisms for the creation of the newly-found "overionized" plasma and MM SNRs in 
general. We conclude that both the plasma structure and ionization conditions can be 
explained by thermal conduction. We finally present the implications of our findings for the 
overall picture of SNR evolution. 

2. Target Selection 

In order to investigate the ionization conditions in MM SNRs, we adopted the following 
method to supplement global NEI model fits. First, we estimate the electron temperature, 
kT e , from the continuum shape of the X-ray spectrum. Next, we measure the intensity ratio 
of the H-like Ka (hereafter Lya) line to the He-like Ka (hereafter Hea) line for the heavy 
elements, which provides a measure of the degree of ionization as the ionization temperature, 
kT z . Comparison between kT e and kT z directly gives the ionization state of the SNR from its 
X-ray spectrum. This method is to produce robust results, then we need to utilize data from 
an X-ray spectrometer that resolves the Lya line from the Hea line. One such instrument 
is the Solid-state Imaging Spectrometer (SIS) onboard ASCA, which has energy resolution 
of 5% (FWHM) at 1.5 keV. However, the SIS energy resolution degraded over the course 
of the mission as a consequence of accumulated radiation damage (Yamashita et al. 1999). 
Therefore, only data from early ASCA observations (from the launch in 1993 through 1994) 
are suitable for this analysis. 

Table 1 tabulates the ASCA observations of candidate MM SNRs from Rho & Petre 
(1998). Among the candidates, MSH 11-54 (G292.2+1.8), their group "F" remnant, is 
excluded because it is considered to be a reverse- shocked, ejecta-dominated young remnant, 
showing strong O and Ne lines with a large velocity of 2,000 km s -1 (Goss et al. 1979; Murdin 
& Clark 1979), and a faint shell in X-rays (Park et al. 2002). For the systematic study of 
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ionization states, we used only the ASCA data. Thus, our criteria of selecting targets in this 
paper are: (1) MM SNR candidates assigned as A to E in Table 1, and (2) SNRs observed 
with ASCA from 1993 to 1994. There are eight SNRs satisfying these criteria. However, the 
data from 3C400.2 and MSH 11-61A, observed in 1994, are not suitable for our analysis as 
mentioned in the notes of table 1. Therefore, we analyzed the following six remnants: IC 443 
(G189. 1+3.0), W49B (G43.3-0.2), W28 (G6.4-0.1), W44 (G34.7-0.4), 3C391 (G31.9+0.0), 
and Kes 27 (G327.4+0.4). 



3. Data Analysis and Results 

For the series of analysis, we extracted all the data from ASCA public archive, and 
screened them using the the NASA/GSFC revision 2 standard data processing criteria. We 
applied no Residual Dark Distribution (RDD) correction because the early-mission SIS data 
show few or no radiation damage effects. As most of these remnants are located near the 
Galactic plane where the Galactic Ridge X-ray emission (GRXE) is not negligible, we used 
source-free regions of the individual observations for background subtraction. Exceptions to 
this background subtraction approach were IC 443 and W49B, for which we used blank-sky 
data: IC 443 is located in the Galactic anti-center region, and the results for W49B are not 
sensitive to the background used (Hwang, Petre, & Hughes 2000). 

When investigating the ionization state, we used a model consisting of a CIE plasma 
(VMEKAL in XSPEC, Mewe, Gronenschild, & van den Oord (1985); Kaastra (1992); Liedahl, 
Osterheld, & Goldstein (1995)) and several Gaussian components. A CIE plasma model 
with metal abundances set to zero represents the continuum shape (i.e., kT e ), while Gaussian 
components are used to measure the line intensity ratios (Lya/Hea) of heavy elements (i.e., 
kT z ). We fixed the Gaussian widths at 10 eV for Hea of Mg, Si, S, and Ar, and 20 eV 
for Hea of Ca, in order to account for the fact that the Hea lines are triplets. The widths 
of all other lines were fixed at 0.1 eV, corresponding functionally to monochromatic lines 
compared with the SIS energy resolution. 

For remnants whose angular size is larger than the SIS field of view (FOV), we extracted 
spectra from a few regions and investigated the plasma structure and ionization state of each. 
Below, we present the results for individual remnants. 
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3.1. IC 443 

The details of the analysis and results for IC 443 are described in Kawasaki et al. (2002); 
we briefly summarize the results here. IC 443 has a ~45' diameter, and we extracted the 
spectra from two regions: a central region where the X-ray surface brightness is highest, and 
a region near the northern rim. Both spectra require two plasma components (1.0 keV and 
0.2 keV); the emission of the 0.2 keV plasma is stronger in the region near the shell than 
the center. In addition, the mean electron densities of the 1.0 keV and 0.2 keV plasmas 
estimated from the emission measures are 1 cm -3 and 4 cm -3 respectively. These results can 
be accounted for by a simple model in which IC 443 has a hot (1.0 keV) interior surrounded by 
a cool (0.2 keV) and denser shell. In the high-temperature plasma, the ionization temperature 
from the Lya/Hea ratio of S is ~1.5 keV, significantly higher than the electron temperature 
of 1.0 keV. Neither an additional, hotter plasma component nor a multi-temperature plasma 
successfully accounts for this feature, and we conclude that the interior hot plasma of IC 443 
is overionized. 

We attempted to verify this result using XMM-Newton data. Unfortunately, the obser- 
vation of IC 443 suffers from high background; the available data do not provide results as 
statistically accurate as those from ASCA. 

3.2. W49B 

W49B appears centrally brightened in X-rays and shell-like at radio wavelengths, like 
other MM SNRs. However, this remnant is thought to be ejecta-enriched, and thus relatively 
young, based on the high (3-5 times as solar) metal abundances (Fujimoto et al. 1995; 
Hwang, Petre, & Hughes 2000). Since W49B is compact, with a size of only 4'x3' we used 
the Chandra images with 0.5" spatial resolution to ascertain the energy dependence of its 
morphology. Figure 1 shows the Chandra ACIS images in the 1.0-2.12 keV, 2.12-4.10 keV, 
and 4.25-6.30 keV energy bands. They show that the highest energy image is more centrally 
concentrated than those in the lower energy bands. The SIS spectrum of the entire remnant 
in the 0.8-10.0 keV energy band (2CCD mode data) requires at least a two-temperature 
(0.24+°;^ keV and 1.701°;°^ keV) CIE plasma (VAPEC in XSPEC, Smith et al. (2001)) with a 
large absorption column density of 5.23^q;Jq x 10 22 cm -3 . Furthermore, we had to add four 
Gaussian components to reproduce unusually strong emission lines of Ar Lja and Ca Lya, 
whose center energies were fixed at theoretical values of 3.324 keV and 4.105 keV respectively, 
as well as Cr Hea and Mn Hea lines (Hwang, Petre, & Hughes 2000). The reduction of x 2 by 
adding each Lya line is more than 10 for 147 d.o.f, and an F-statistic test gives a probability 
greater than 97% indicating that these lines are real. The resultant % 2 /d.o.f. is 233/146 (see 
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Tab. 2 and Fig. 2). 

We therefore evaluated the ionization temperatures of Ar and Ca and compared them 
to the continuum temperature. In order to achieve a better statistical accuracy, we added 
the two data sets taken in 1CCD mode (13ksec and 19ksec SIS data) and fitted them simul- 
taneously. We used the spectra in the 2.75-6.0 keV energy band since the contribution of the 
low-temperature plasma is negligible at these energies. We applied a model consisting of the 
CIE plasma (VMEKAL) and 9 Gaussian components (Hect, Lya, He/5, He7 of Ar, Hea, Lya, 
He/3, He7 of Ca, and Hea of Cr) as shown in Figure 3 (left). Line center energies of Gaussian 
components were fixed at the theoretical values except for Ar Hea, Ar Lya, Ca Hea, and 
Ca Lja. For the W49B plasma with strong emission lines, line blending can be a significant 
issue when measuring line fluxes with a moderate resolution spectrometer like the SIS. A 
non-negligible flux from S H-like K/3 (Ly/3) at 3.107 keV can blend with Ar Hea (~3.13 keV), 
as can He/5 and He7 emission of S at 2.88 keV and 3.02 keV respectively. In these fits, we 
fixed the S abundance to 3.9 solar as derived from the 0.8-10.0 keV spectrum fitting, instead 
of adding further Gaussian components. We ignored the Ar Ly/5 (3.936 keV) flux, which 
could mildly contaminate (5-10%) the Ca Hea. For He/5 and He7 line blending, we modeled 
the He7 line with an intensity relative to the He/5 line fixed at 0.6 (Hwang, Petre, & Hughes 
2000). Both Ar and Ca, the Lya/Hea ratio requires an ionization temperature of ~2.5 keV; 
this is higher than the ~1.8 keV electron temperature (see Fig. 3 (right) for the confidence 
contours). Adding a third plasma component with higher (> 1.8 keV) temperature can in 
principle account for the strong H-like lines of Ar and Ca, but it is required to have an im- 
plausible combination of small emission measure and extremely large (> 100) abundances. 
Therefore, we suggest that the high-temperature plasma of W49B is "overionized" , similar 
to IC 443. W49B is thus the second SNR showing overionization. 

Our results may be compared with the detailed spectral analysis by Hwang, Petre, & 
Hughes (2000). They showed the constraints on kT z for the Lya/Hea of Ar and Ca to be 
2.2-2.7 keV. The values are higher than the kT c of 1.7 keV inferred from the bremsstrahlung 
continuum. Our results are consistent with theirs, only our interpretation differs. 

Taken together, the Chandra images and ASCA spectra indicate the presence of central 
hot and outer cool plasmas. We thus suggest that W49B consists of two typical plasma 
components (a 1.7 keV hot interior and a 0.24 keV cool shell), like that of the Sedov model, 
though Fujimoto et al. (1995) suggested a possible stratification of ejecta. 

The electron densities can be estimated from the emission measures to be 
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for the high- and low-temperature components respectively, assuming that the two compo- 
nents are in pressure equilibrium and the X-ray emitting volume extends to the radio shell 
(which we approximate as a spherical region with angular radius 9=2' ). Here, f\ and f 2 
(fi + f2 — 1) are the filling factors of these components and d is the distance to W49B. This 
plasma structure, temperature and density distributions, is consistent with the idealization 
of the Sedov model. 

The X-ray emitting masses of the components can also be estimated: 
Mi 

M 2 

The masses of the high- and low-temperature plasmas are much higher than that of any 
plausible progenitor, even if d = 8 kpc which is favored by Moffett & Reynolds (1994). 
We should note here that the masses are estimated assuming uniform density distributions 
within the remnant. The Chandra images show some clumpy structures as well as the ex- 
tended components. In fact, the mass of the high-temperature plasma could be ~ M if the 
X-rays come only from clumps (Hwang, Petre, & Hughes 2000). However, the extent of the 
1.00-2.12 keV Chandra image (see Fig. 1) indicates that the mass of the low-temperature 
component is much higher than the ejecta mass. Coupled with the surprisingly close proxim- 
ity of the plasma to ionization equilibrium (in stark contrast to ejecta-dominated remnants 
(e.g.,G292.0+1.8 - (Hughes & Singh 1994)), the high mass leads us to suggest that W49B 
is not a young, ejecta-dominated, but an evolved MM SNR. 

The mass and ionization state of W49B are reasonably explained if the plasma is dom- 
inated by the shocked ambient medium, but its large abundances and high temperature 
(up to ~1.7 keV) strongly suggests the the presence of ejecta. Overabundances might be ex- 
plained if W49B exploded in the massive (~10 2 M ) circumstellar cloud with large amount of 
heavy elements seeded by previous SNRs. The X-ray measured absorption column density is 
5.23^o}qX10 22 cm -3 , significantly higher than the Galactic column density of 4.4xl0 22 cm -3 
measured in this direction using the 21cm and CO observations (Dickey & Lockman 1990; 
Stark et al. 1985). The X-ray absorption is mainly due to heavy elements, not hydrogen. 
Thus, this high column density might indicate that W49B is embedded in an unusually metal 
rich environment. More detailed investigation using higher resolution spectrometers will be 
needed to resolve this question, and it is clearly outside the scope of this work. 
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3.3. W44 

We show in Figure 4 the smoothed images of W44 in the 0.7-1.5 keV and 1.5-4.0 keV 
bands using the Gas Imaging Spectrometer (GIS). The soft (0.7-1.5 keV) image has a larger 
angular extent than the hard (1.5-4.0 keV) one, suggesting spectral variations due to a 
temperature gradient. Furthermore, clear differences are found in fits to spectra from the 
central and northern regions shown in Figure 4. The temperature of the central region, 
0.84+q;q5 keV, is higher than that of the northern, 0.64l[j;o5 keV. This is consistent with 
the fact that the hard band image is more centrally concentrated. We also found higher 
abundances of Mg and Si, significant at the 90% confidence level, in the central region: 
e.g., the abundance of Si is 1.18lo' 12 m the cen ter while in the north it is 0.50TJoo4- This 
enhancement could be due to residual ejecta, as first detected in the Cygnus Loop (Miyata 
et al. 1998). The best fit ionization timescale is >10 12 s cm -3 for both regions. Details 
are tabulated in Table 3. The Chandra data also show the remnant's hot and metal-rich 
projected center with bright knots (Shelton, Kuntz, & Petre 2004a). Note that the SIS data 
exhibited a slight energy shift, which are also shown in the Chandra data (Shelton, Kuntz, 
k Petre 2004a). 

The ionization state of the W44 plasma was investigated using line intensity ratios. In 
this fit, we used the 1.6-6.0 keV band spectra (see Fig. 5(a)). An interstellar absorption 
value was fixed to 8.9 x 10 21 cm~ 2 , which was obtained from the NEI fits. The ionization 
temperature of Si in the central region, obtained from the fitted Lya/Heo; ratio, is ~ 0.7 keV. 
This is consistent with the electron temperature of ~ 0.8 keV as shown in Figure 5(b), as 
well as with the NEI model fitting. The inferred kT z of Si in the north region is ~ 0.8 keV, 
with a large (0.5-1.5 keV) statistical uncertainty. 

3.4. W28 

The X-ray morphology of W28 is basically center-filled, but is more complex than the 
other MM SNRs. It has ear-like segments of limb-brightened shell in the northeast and 
northwest (Rho & Borkowski 2002). In addition, we found with the ASCA GIS that the 
emission from the northeastern part is strong in the soft band (0.5-1.5 keV) but not in the 
hard (1.5-4 keV). This is also indicated by Figures 6 and 8 in Rho & Borkowski (2002). 

We extracted the spectra from the central area, named "center" (a rectangular region 
with 8'x7' centered at R.A.=18 h 00 m 26!, Dec.=-23 26'49") and the eastern area, named 
"east" (an elliptical region with 4.8'x4.1' centered at R.A. = 18 h 00 m 36!, Dec.=-23°23'45"). 
The spectrum was well fitted by the two NEI (VNEI) components with temperatures of 
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0.6 keV and 1.4 keV in each region (see Tab. 4). The ionization timescale of both components 
are around 10 12 s cm" 3 . The abundances of Ne, Mg, Si, S, which were linked together for the 
two plasmas, are fitted with <1 times solar values (Anders & Grevesse 1989). The column 
density is 5.1xl0 21 cm" 2 and 5.4xl0 21 cm" 2 for the center and the east, respectively. We 
found a larger fraction of high-temperature plasma in the center region (the emission measure 
ratio is 0.5llo; 2 g) than that in the east (the ratio is 0.3llo!o4)- m contrast, the northeast- 
ern shell consists of only a low-temperature plasma (A;T C =0.56 keV, n c t=1.7xl0 13 s cm -3 ) 
(Rho & Borkowski 2002). Therefore, we suggest that low-temperature component becomes 
dominant toward the northeastern rim of the remnant where dense molecular clouds are 
present. These results are consistent with those of IC 443 (Kawasaki et al. 2002), implying 
that the plasma structure of a hot interior and a cool exterior is present in W28. On the 
other hand, the southwestern shell region shows hard X-ray emission. Rho & Borkowski 
(2002) showed that the spectrum in this region can be fitted with a temperature of 1.5 keV 
and n e t=1.5xl0 11 s cm" 3 . This implies that the plasma structure of W28 is more complex 
than allowed by the above simple model. 

The ionization state of the high-temperature plasma was examined using line intensity 
ratios. We used the spectra in the 1.65-6.0 keV energy band although they include con- 
tamination by the low-temperature plasma (see Fig. 5(c)). We thus modeled the column 
density and low-temperature plasma with all parameters fixed to the values estimated from 
the component fits. The ionization temperature of Si in the east region, from the fitted 
Lya/Hea ratio, is close to the fitted electron temperature as shown in Figure 5(d), although 
the poor statistics do not constrain the temperatures well. The inferred kT e and kT z of Si 
in the central region are 1.19^Q; 2 g keV and >1.06 keV, respectively. These results indicate 
that the high-temperature plasma is near CIE in each region. Hence, our direct estimate 
of the ionization state suggests that the plasma in W28 has reached ionization equilibrium, 
consistent with the results of the NEI fits. 



3.5. 3C391 

3C391 shows an elliptical X-ray morphology, with a size of 7'x5', and major axis oriented 
from northwest to southeast. Chen & Slane (2001) found that the ASCA hard band (2.6- 
10 keV) image shows stronger emission in the northwest, in anti-correlation with the soft band 
(0.5-2.6 keV) image. This difference is due to the larger column density in the northwest. 
The variation of the column density across the remnant is in agreement with the presence of 
a molecular cloud to the northwest. No temperature gradient has been detected across the 
remnant, a result confirmed by Chandra (Chen et al. 2004). 
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We extracted the spectrum of the entire remnant, in the energy band 0.8-5.0 keV. 
We first applied a single NE1 model, and obtained an electron temperature and ionization 
timescale 0.53±g;| keV and 2.5 (>0.9) xlO 12 s cm . These suggest that the plasma has 
reached ionization equilibrium (see Tab. 5). We then measured the ionization degree of 
the plasma using line intensity ratios (see Fig. 5(e)). Since the spectrum shows prominent 
emission from Hea lines of Mg, Si, and S with good energy resolution, we left their center 
energies free; those of the others were fixed to the theoretical values. Figure 5(f) shows 
confidence contours of kT z to kT e for Mg and Si, and indicates that the ionization temperature 
is comparable with the electron temperature. Thus, both our direct estimate and NEI model 
fitting suggest that the plasma in 3C 391 is near ionization equilibrium. 



3.6. Kes 27 

Kes 27 shows diffuse, thermal X-ray emission with several unresolved point-like sources 
(Seward, Kearns, & Rhode 1996). Using ASCA data, Enoguchi et al. (2002) found a temper- 
ature gradient from the rim (0.59l[};{£ keV) to the interior (0.84±0.08 keV). This temperature 
variation is similar to that in W44 (see §3.3). 

The 0.8-6.0 keV SIS spectrum of the entire remnant is well fitted by a single NEI model 
(see Tab. 6). The inferred electron temperature, column density, and ionization timescale 
are 0.9±0.1 keV, 2.2±0.1xl0 22 cm" 2 , and > 3 x 10 11 s cm" 3 , respectively. We measured the 
ionization state of Kes 27 using the line intensity ratios of Si and S, as shown in Figure 5(g) 
and (h). The center energies of all the lines are fixed to the theoretical values for a 0.9 keV 
CIE plasma. The inferred Lya/Hea ratio for Si indicates that the ionization temperature 
of Si is close to the electron temperature. However, that of S may require an ionization 
temperature larger than the continuum temperature, suggestive of an overionized plasma, 
although this statement is significant at a confidence level of <90%. The results indicate that 
the plasma in Kes 27 is basically in ionization equilibrium, but the strong S Lja emission 
suggests possible overionization. 
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4. Discussion 
4.1. Comparison with Shell-like SNRs 

4-1.1. Plasma Structure 

We first summarize in Table 7 the spatial distributions of electron temperature and den- 
sity for the six MM SNRs. IC 443 and W49B show clear evidence of temperature and density 
gradients; they obviously indicate that these SNRs have a hot interior surrounded by a dense, 
cooler exterior plasma, qualitatively consistent with that of a standard shell-like SNR. We 
note that these temperature and density gradients are not as large as those predicted by 
the Sedov solution. Additionally, W44, W28, and Kes 27, show temperature gradients; they 
have higher electron temperature towards the center where the X-ray surface brightness is 
highest than towards the edge. A similar temperature gradient has been detected from an 
ISM dominated shell-like SNR, the Cygnus Loop (Miyata et al. 1994, 1998), as tabulated in 
Table 7. These results suggest that MM SNRs have experienced similar evolution to shell- 
like SNRs. 3C391 (Chen & Slane 2001; Chen et al. 2004) is the only SNR with no detectable 
temperature or density gradient of the six studied here. 

4-1.2. Ionization States 

We summarize in Table 8 the electron and the ionization temperatures of the MM SNRs 
we have analyzed. Two of them, IC 443 and W49B, show significantly higher ionization 
temperatures than electron temperatures, suggesting the plasma is overionized. For the 
other remnants, the two temperatures are comparable, suggesting that the plasma has almost 
reached CIE. This is supported by the long ionization timescales obtained by the NEI model 
fits. 

For comparison, we tabulate in Table 9 the electron temperature (kT c ), ionization 
timescale (n e t), and ionization temperature (kT z ) of prominent shell-like SNRs. Cassiopeia A 
(Cas A), Kepler's SNR, and Tycho's SNR exemplify the young (~ 500 yr) SNRs, and the 
others are examples of middle-aged (> 1000 yr) SNRs. We inferred kT z of each SNR by the 
following method: from kT c and n e t, the emissivities of Hea and Lya lines for heavy elements 
are calculated using an NEI plasma code (Masai 1984, 1994); their ratio, Lya/Hea, yields 
the ionization temperature kT z . We estimated the Lya/Hea ratios for Ne, Mg, Si, and S 
using this method; the inferred kT z values are consistent with one another within ~0.2 keV. 
The validity of this estimate is checked with the Cas A spectrum as follows. We fitted the 
SIS spectrum of Cas A with two NEI (VNEI) models; the measured kT e and log(n c t) of the 
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high-temperature plasma are 3.02 keV and 10.83, respectively, from which we estimate kT z 
of S as 1.071^4 keV. On the other hand, we derived the kT z of S for the high-temperature 
plasma to be 0.94 ±0.07 keV, from the fitted intensity ratio of Lja to Hea (Lya/Hea) as we 
did for MM SNRs. In this fit, the low-temperature plasma was modeled as another plasma 
component, in the same way as for W28. The ionization temperatures derived with the two 
methods are consistent within 0.1-0.2 keV. Thus, we confirmed that the values tabulated in 
Table 9 properly represent kT z of these shell-like SNRs, and can be compared with those for 
MM SNRs. 

Figure 6 shows the correlation between kT c and kT z for both the MM and shell-like 
SNRs. We can see a systematic difference between the two types of SNRs. All MM SNRs 
show kT z > kT c , which means that the plasma has reached CIE or a state of overionization. 
On the other hand, all of the shell-like SNRs show kT z < kT e , indicating underionized plasma 
conditions. These results strongly suggest that MM SNRs are systematically more ionized 
than shell-like SNRs. 



4.2. Formation Mechanisms of the Overionized Plasma 

We also found that Lja line emission from IC 443 and W49B is unusually strong com- 
pared with Hea, indicating that the interior plasma in these remnants is overionized. In the 
low density plasma found in SNRs, we are accustomed to encountering underionized plasma, 
like that detected in shell-like SNRs, because the ionization timescale is typically long com- 
pared with the time since the bulk of the gas was shock-heated. Nevertheless, conditions 
can exist in which such SNR plasmas can become overionized. 

We consider two possible causes for overionization: photoionization, and rapid gas cool- 
ing compared with recombination. For photoionization, a low-density SNR plasma itself 
does not produce sufficient radiation to photoionize heavy elements, and there is no strong 
X-ray source in any of these remnants (see Fig 1 for W49B, and Kawasaki et al. (2002) for 
IC 443.) On the other hand, there are three possible cooling mechanisms in a gas: radia- 
tion, expansion, and thermal conduction. Kawasaki et al. (2002) compared these timescales 
with that of recombination, and showed that in IC 443 thermal conduction could drive the 
plasma into an overionized state. A similar comparison can be carried out for W49B. The 
conduction timescale, 9 x 10 10 (n 1 /2.5 cm- 3 )(/ T /10 19 cm^ifeTi/l.S keV)- 5 / 2 (ln A/32.2) sec, 
is shorter than the recombination timescale of 4 x 10 11 (ni/2.5 cm -3 ) sec, while the radia- 
tive cooling timescale is 3 x 10 14 (A;Ti/1.8 keV)(#i/1.7') 3/2 (d/11.4 kpc) 1 / 2 sec. Here, n u U 
and kT\ are the number density, the radius in arcminutes, and the temperature of the inner 
plasma of W49B, respectively, and It, and In A are the scale length of the temperature gradi- 
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ent and Coulomb logarithm. We note that the cooling timescale via expansion is 12.5 times 
the remnant's age (Masai 1994), meaning that the remnant has not been cooled enough at 
present time. This indicates that thermal conduction is the dominant process producing the 
overionized plasma in W49B as well as in IC 443. 

The presence of a magnetic field will reduce the heat flux in principle, because diffusion 
perpendicular to the field direction is suppressed. The magnetic field reduces the heat 
flux by a factor of ~ (cos 2 6) = 1/3 if it is tangled, but the conduction timescale still 
remains comparable to the recombination timescale. Although detailed information about 
the magnetic field in W49B is unavailable, we conclude from the above argument that it is 
possible for the gas to be overionized due to thermal conduction. The fact that we have not 
detected in any SNR as large a temperature distribution as predicted by the Sedov solution 
further suggests that thermal conduction is operating in many remnants. 

The process of thermal conduction should arise in all SNRs, as Cox et al. (1999); Shel- 
ton et al. (1999) have introduced and discussed thermal conduction over the lifetime of an 
SNR. Therefore, we investigated the differences between overionized SNRs and the other 
MM SNRs by comparing the ionization states to the parameters related to thermal con- 
duction. Figure 7 (left) show the ionization states (i.e., kT z /kT e ) plotted against the ratios 
between the conduction timescale and the recombination timescale (icond/Womb) for the six 
remnants we have analyzed. In the calculation of the conduction timescale, we ignored the 
effect of the magnetic field. The ratio of the timescales has strong dependence on the density, 
temperature gradient, and electron temperature as t C ondArccomb oc n 2 l^(kT c )~ 5 ^ 2 . This figure 
suggests that conduction in W44 and 3C391 is not strong enough to produce overionization, 
consistent with the observations. On the other hand, our estimates suggest Kes 27 and W28 
might be overionized, even though our analysis reveals no clear evidence. Two possibilities 
are considered for the discrepancy: one is that we could not detect the evidence because 
of the poor statistics as described in §3.6 for Kes 27; the other is that the magnetic field 
or some other process reduces the heat flux. We also compared the ionization state to the 
temperature gradient within the remnants because the temperature gradient strongly affects 
the thermal conduction as described above. In Figure 7 (right) we compare the ionization 
state against the temperature ratio between the exterior and the interior (kT e (exterior) /kT e 
(interior)) for the six MM SNRs. This figure indicates that overionized SNRs have a larger 
temperature gradient than the other MM SNRs(Note, however that temperatures of the in- 
terior and the exterior were not estimated in the same way for the six remnants.) We also 
suggest that MM SNRs with no overionization, such as W44 and 3C391, can be those which 
conduction has taken further from an overionized state since thermal conduction gradually 
reduces the temperature gradient within the remnant. 
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4.3. SNR Evolution and Formation of MM SNR 

As we described in §4.2, thermal conduction arises in an SNR, and it equilibrates the 
temperatures while keeping the pressure balanced. It suggests that the hot interior plasma 
gradually becomes denser. Therefore, we propose the following unified scenario for SNR 
evolution. We present in Figure 8 a schematic view on the evolution of the ionization state 
in a remnant. Numerical simulations by Shelton (1998, 1999) also predict both the center- 
filled X-ray morphology and overionization as a result of the SNR evolution. 

In a young remnant at an early adiabatic or Sedov phase, the strong front shock heats 
the ambient ISM to high temperatures and bright X-rays appear along the outer rim. The 
hot gas within the post-shock interior is rarefied compared to the shell regions, and hence, the 
"shell-like" morphology in X-rays is formed. The observed plasma is mainly from recently 
shocked regions where ions are being ionized to higher ionization states by electrons on the 
timescale of n e t ~ 10 12 s cm -3 . Therefore, an underionized plasma is detected (phase- 1 in 
Figure 8). 

As the remnant evolves, the hot interior plasma gradually approaches a uniform tem- 
perature via thermal conduction, while maintaining the pressure equilibrium. Thus, the 
interior also becomes denser, resulting in an enhancement of the X-ray emission from the 
central region. In contrast, the shock speed decreases and the X-ray emission from the shock- 
heated shell becomes too soft to pass through the ISM. These lead the X-ray morphology to 
change from "shell-like" to "center-filled" on the conduction timescale, which is the forma- 
tion mechanism of MM SNR. Since the conduction timescale is roughly comparable to the 
ionization timescale, t con d ~ trecomb — 10 12 s cm -3 , ionization states of the remnants reach 
nearly an ionization equilibrium as derived from our analysis of MM SNRs. This phase can 
be termed the "conduction phase" since thermal conduction plays an important role to the 
SNR plasma. 

If the conduction timescale is shorter than the ionization timescale, the interior plasma 
could become "overionized" over a certain interval, as observed in IC 443 and W49B (phase-2 
in Figure 8). It is uncertain how long this interval lasts, but the plasma will eventually reach 
another ionization equilibrium because conduction becomes insufficient as the interior plasma 
cools and the temperature gradient disappears (phase- 3 in Figure 8). When the conduction 
becomes fully efficient, no temperature gradient remains in the remnant, as observed in 
3C391. Figure 7 (right) may indicate the evolution from the overionized SNR (phase-2) to 
the fully conducted SNR (phase-3). 

This interpretation is supported by the X-ray morphological analysis of SNRs in the 
Magellanic Clouds that "centrally brightened" SNRs are larger in size than "shell-type" SNRs 
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(Williams et al. 1999), and by G65.3+5.7, the remnant newly identified to be in transition 
from shell-like to MM (Shelton, Kuntz, & Petre 2004b). It is controversial to include W49B 
to the MM SNR type since it has been considered to be ejecta-dominated plasma. However, 
from the point of view of ionization, thermal conduction plays an important role in W49B, 
whether it is ejecta or not, and it supports the above interpretation. 

5. Conclusions 

We have discovered the following results from the systematic analysis of the six MM SNRs. 

1. Hot interior and cool outer plasma components are detected in five MM remnants. A 
higher density is found in the cool exterior than that in the hot interior from both 
IC 443 and W49B. These results suggest a plasma structure consisting of a central 
hot region surrounded by a cooler, denser shell, similar to the expected structure of a 
shell-like SNR. 

2. The X-ray emitting gas in IC 443 and W49B is "overionized;" the gas in the other 
MM SNRs are found to be near ionization equilibrium. As a consequence, the degree 
of ionization in MM SNRs is systematically higher than in shell-like SNRs. 

The first result indicates that MM SNRs have experienced a similar evolution to shell-like 
SNRs, and the second result indicates that MM SNRs are more evolved than shell-like SNRs. 
We suggest that thermal conduction plays an important role in producing the overionized 
interior plasma in IC 443 and W49B. It also acts to diminish the temperature and den- 
sity gradients within all remnants, leading naturally to a "center-filled" X-ray morphology. 
Hence, we conclude that MM SNRs are those which have evolved over a conduction timescale 
of ~ 10 11 " 12 sec with thermal conduction active within the remnants. It is also suggested 
from the fact that plasmas of MM SNRs are nearly in CIE, with the estimated ionization 
timescale of ~ 10 12 sec. We propose to call this phase the "conduction phase" since thermal 
conduction drastically changes the plasma structure of an SNR in this phase far from that 
in the Sedov phase. By introducing thermal conduction, we showed that MM SNRs are not 
special SNRs but those in a certain stage of an SNR evolution. This leads to the construction 
of a unified model of SNR evolution. 

We thank Dr. Manabu Ishida and Dr. Kuniaki Masai for helpful discussion, and Dr. 
Kazuo Makishima for thoughtful comments and suggestions. 
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Fig. 1.— Chandra ACIS images of W49B in (left) 1.00-2.12 keV, (middle) 2.12-4.10 keV, 
and (right) 4.25-6.30 keV bands. Grayscales are linear and the contour levels in each figure 
are 20%, 40%, 60%, and 80% of the peak surface brightness. Offset center in each figure is 
R.A.=19^11P1!35, Dec.=9°06 / 7.64 // (J2000). All images are derived from Chandra Supernova 
Remnant Catalog Web page (http://hea-www.harvard.edu/ChandraSNR/). 
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Fig. 2. — ASCA SIS spectra of the entire region of W49B with the best-fit 2 CIE plasma 
and Gaussians model. Dashed, dash-dotted, and dotted lines represent the low-T plasma, 
high-T plasma, and additional Gaussian components respectively. The lower panels show 
the residuals of the fits. 




Fig. 3.— left ASCA SIS spectra of W49B in 2.75-6.0 keV energy band with the best-fit 
model of VMEKAL (dashed) and 9 Gaussian lines (dotted). The data in 2CCD mode are 
plotted in black, and those in 1CCD mode are in red and blue. The column density is fixed 
to ]Vh=5.23x 10 22 cm -2 . The minimum value of x 2 /d.o.f. is 210/235. The lower panel shows 
the fit residuals. There are residuals at ~3.5 keV, which can be identified as a He-like Ka 
line of potassium, right: Confidence contours of ratio of the ionization temperature (kT z ) to 
the electron temperature (kT e ) for Ar (black) and Ca (red) with a fixed column density of 
5.23 xlO 22 cm -2 . The confidence levels are 99%, 90% and 67%. The dashed line represents 
CIE. 
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Fig. 4.— Smoothed GIS images of W44 in the 0.7-1.5 keV and (left) 1.5-4.0 keV (right) 
bands. Grayscales are in linear scale and contour levels in each figure are 20%, 40%, 60%, 
and 80% of the peak surface brightness. Spectral extraction regions named "center" and 
"north" are shown in white ellipse and circle, respectively. 
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(a) W44 spectrum 1.6-6.0 keV 



(b) W44 contours: center region 
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Fig. 5. — X-ray spectra and confidence contours of the ratio of the ionization temperature 
(kT z ) to the electron temperature (kT e ) for W44, W28, 3C391, and Kes 27. In each spectrum, 
the dashed line (and dash-dotted line for W28) shows the continuum component while dotted 
lines represent line features. The lower panel shows the residuals of the fit. In each contours, 
confidence levels are 99%, 90% and 67%. The dashed line represents kT e = kT z , meaning 
that the plasma is at the CIE state. 
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Fig. 6. — The ionization temperature vs. the electron temperature for various MM SNRs (red 
open square) and Shell-like SNRs (blue open circle). The dashed- line represents kT e =kT z , 
implying that the plasma is at the CIE state. 
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Fig. 7. — Comparisons between overionized SNRs and the other MM SNRs. left: The rela- 
tions of the ionization states (kT z /kT e ) to the timescale ratio (t C ondArecomb) for the MM SNRs. 
The typical error of the ratio due to uncertainties in It is also shown at the lower left of the 
panel. The condition t C ondArccomb > 1 means that it is hard for the plasma to be overionized 
due to thermal conduction, right: The relations of the ionization states (kT z /kT e ) to the 
electron temperature gradients (kT e (exterior) /kT e (interior)) for the MM SNRs. The dashed 
line denotes kT e (exterior) /kT c (interior) = 1, meaning a uniform temperature within the 
remnant. For 3C391 with no detection of temperature difference within the remnant, we 
assumed the uniform temperature. 
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Fig. 8. — Schematic view of the evolution of the ionization state in an SNR along time. The 
horizontal line indicates the time normalized by the plasma density, and the vertical line 
qualitatively shows the electron and ionization temperatures. The ionization temperature 
follows the electron temperature with a delay of n e t ~ 10 12 s cm~ 3 . 
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Table 1. The list of MM SNR candidates with ASCA observation logs 



SNR name 


Size (arcmin) 


Group a 


ASCA Obs. 


Exposure time b 


W28 


42 


A 


1994/4, 1995/9 


50ks(1994/4, 4CCD) 


W44 


28x34 


A 


1994/4, 10 


24ks(1994/4), 26ks(1994/10) 


3C400.2 


28 


A 


1994/10 c , 1996/4 




Kes 27 


20 


A 


1994/8 


23ks 


MSH 11-61A 


12 


A 


1994/3 d , 1995/1 




3C391 


8 


A 


1994/4 


44ks 


CTB 1 


34 


A 


1996/1 




W51C 


25x19 


B 


1995/10 




CTA l c 


80 


B 


1996/1, 1997/2 




VV \JO 




R 

1 J 






HB21 


120x90 


B 


1997/6 




IC 443 


50 


C 


1993/4, 1994/3, 1998/3, 1999/3 


72ks(1993/4) 


Kes 79 


11 


C 


1995/4 




HB3 


70 


c 


1996/8 




G327.1-1.P 


14 


D 


1996/3 




CTB 104A 


80 


D 


no observation 




W49B 


3.5 


E 


1993/4, 5, 10, 11 


46ks(1993/4) 


3C 397 


6 


E 


1995/4 




MSH 11-54 


10 


F 


1993/8 





a classification of the MM SNR candidates assigned by Rho & Petre (1998): The group A 
remnants are the prototypical MM SNRs, group B remnants are probably MM SNRs, and 
group C remnants show additional spatial components. Group D remnants are those whose 
data quality are too poor to allow any determination. Group E are probably not MM SNRs; 
nevertheless, they have similar characteristics, and so their classification remains uncertain. 
Group F is definetely not an MM SNR. 

b Effective SIS (SIS0 + SIS1) exposure time after screening are listed for the targets that 
we analyzed in this paper. 

c The target was out of the field of view because of a pointing error 
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d Due to observation at small Earth angle, only a small fraction of the SIS data remains 
after the screening (Slane et al. 2002). 

C ASCA observations revealed the non-thermal component that suggests these remnants 
are a true composite type (Slane et al. 1997; Sun, Wang, & Chen 1999). 
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Table 2. The best-fit parameters of the two CIE model for the W49B spectra 





parameter (unit) 


value 


(90% cr.) 


wabs a 


N u (10 22 cm" 2 ) 


5.23 


(5.13-5.34) 


Low-T plasma 


kT e (keV) 


0.24 


(0.22-0.28) 




E.M. b 


2.7 


(1.0-5.8) 


High-T plasma 


kT e (keV) 


1.70 


(1.66-1.72) 




E.M. b 


0.094 


(0.090-0.097) 




Mg 


0.55 


(0.18-0.92) 




Al 





(<1.74) 




Si 


3.5 


(3.2-3.8) 


Abundance 


s 


3.9 


(3.7-4.3) 




Ar 


3.2 


(2.8-3.7) 




Ca 


4.5 


(4.0-5.1) 




Fe 


5.5 


(4.9-6.2) 




Ni 


16 


(7.4-25) 


Ar XVIII Ka 


flux c 


7.5 


(4.4-10.7) 


Ca XX Kf3 


flux c 


5.2 


(3.4-7.1) 


Cr XXIII Ka 


flux c 


2.0 


(0.7-3.4) 


Mn XXIV Ka 


flux c 


1.0 


(<1.6) 



X 2 /d.o.f. 233/146 



Note. — Each abundance is relative to the solar value and those of two plasma components 
are linked together. The 90% confidence range (cr.) for temperature, abundances, and line 
parameters are given in brackets. The line width of each Gaussian profile is fixed to 0. 

a The column density of an interstellar absorption assuming solar abundances. 

b Emission Measure: the unit is 4ird 2 x 10 14 cm" 5 where d is the distance to W49B. 

c The unit is 10~ 5 photons cm" 2 s _1 . 
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Table 3. The best-fit parameters for the W44 spectra 





Center 


North 


parameter (unit) 


value (90% c.r.) 


value (90% c.r.) 


N u (10 22 cm- 2 ) 


0.89 (0.84-1.04) 


1.03 (0.90-1.25) 


kT c (keV) 


0.84 (0.79-0.86) 


0.64 (0.59-0.68) 


n c t (10 12 cm" 3 s) 


10 (>1.1) 


9.5 (>1.0) 


E.M. a 


6.1 (5.4-7.8) 


10.5 (8.4-13.5) 


Abundance 






Ne 


0.60 (0.16-0.82) 


0.33 (0.16-0.58) 


Mg 


0.78 (0.51-0.89) 


0.21 (0.13-0.30) 


Si 


1.18 (0.96-1.32) 


0.50 (0.46-0.60) 


S 


0.82 (0.61-1.01) 


0.50 (0.33-0.65) 


Fe 


(<0.04) 


0.02 (<0.10) 


Ni 


(<0.2) 


(<0.4) 


xVd.o.f. 


91/87 


83/83 



Note. — Each abundance is relative to the solar value. The 90% confidence range (c.r.) 
for temperature, ionization timescale, and abundances are given in brackets. 

a Emission Measure: the unit is 4ird 2 x 10 12 cm -5 where d is the distance to W44. 
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Table 4. The best-fit parameters for the W28 spectra 





Center 


East 


parameter (unit) 


value (90% c.r.) 


value (90% c.r.) 


iVH (10 cm ) 


5.1 (4.9-5.9) 


5.4 (5.1-6.1) 


T m 1 /tthti—i-t\ 

Low-T plasma (VNEI) 






kT (kev) 


0.62 (0.59-0.65) 


0.61 (0.60-0.63) 


j- ( 1 nl2 —3 \ 

n e t (10 cm J s) 


31 (>2.7) 


15 (>1.0) 


E.M. 


1.7 (1.2-2.2) 


0.9 (0.8-1.1) 


High-T plasma (VNEI) 






kT (keV) 


1.41 (1.31-1.67) 


1.37 (1.31-1.49) 


n c t (10 12 cm -3 s) 


0.5 (0.3-2.1) 


0.4 (0.3-0.8) 


E.M. ratio c 


0.51 (0.32-0.74) 


0.31 (0.27-0.37) 


Abundance 






Ne 


0.28 (<0.88) 


(<0.52) 


Mg 


0.94 (0.65-1.29) 


1.17 (1.04-1.58) 


Si 


0.50 (0.37-0.58) 


0.65 (0.55-0.77) 


S 


0.10 (<0.29) 


0.54 (0.30-0.83) 


Fe 


0.53 (0.48-0.58) 


0.62 (0.56-0.79) 


xVd.o.f. 


137/129 


162/180 



Note. — Each abundance is relative to the solar value and those of two plasma components 
were linked together. The 90% confidence range (c.r.) for temperature, ionization timescale, 
and abundances are given in brackets. 

a The column density of an interstellar absorption assuming solar abundances. 

b Emission Measure: the unit is And 2 x 10 12 cm -5 where d is the distance to W28. 

c This factor represents the ratio of the high-T emission measure to the low-T one. 



Table 5. The best-fit parameters for the 3C 391 spectrum 



Parameter [unit] Entire region (90% c.r.) 



N H [10 22 cm" 2 ] 


2.99 (2.89-3.10) 


kT c [keV] 


0.53 (0.50-0.57) 


n e t [10 12 cm- 3 s] 


2.5 (> 0.9) 


Mg 


1.24 (1.03-1.48) 


Si 


1.12 (0.98-1.26) 


S 


0.81 (0.56-1.07) 


E.M. 11 


14.4 (11.9-17.8) 


X 2 /d.o.f 


93/83 


Note. — Each abundance is relative to the solar value. The 90% confidence range (c.r.) 


are given in brackets. 




^Emission Measure: the unit is And 2 x 10 12 cm 5 where d is the distance to 3C391. 


Table 6. 


The best-fit parameters for the Kes 27 spectrum 


Parameter [unit] Entire region (90% c.r.) 


N H [10 22 cm" 2 ] 


2.2 (2.0-2.3) 


kT c [keV] 


0.94 (0.84-1.01) 


n e t [10 12 cm- 3 s] 


0.7 (> 0.3) 


Mg 


1.0 (0.6-1.8) 


Si 


1.4 (1.3-2.0) 


S 


1.2 (0.8-1.6) 


E.M. 11 


3.09 (2.52-3.40) 


X 2 /d.o.f 


64/64 



Note. — Each abundance is relative to the solar value. The 90% confidence range (c.r.) 
are given in brackets. 

^Emission Measure: the unit is And 2 x 10 12 cm -5 where d is the distance to Kes 27. 
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Table 7. Electron temperatures and densities of mixed- morphology SNRs and Cygnus 

Loop 



SNR 


region 


electron temperature 
[keV] 


density 
[cm" 3 ] 


IC 443 


interior 


1.08 ±0.03 


1.0 




exterior 


0.19 ±0.02 


4.2 


W49B 


interior 


1 70+ - 02 

l - ' u -0.04 


2.5 




exterior 


24 +0 ' 04 


18 


W28 


center, southwest 


1 41+0.25 


0.25« 




center, northeast 


0.62 ±0.03 




W44 


center 
north 


0.841™ 

o-64±8:gg 


0.4( 2 ) t 


3C391 


entire 


C )S +0 ' 04 

U - JO -0.03 


1.1 


Kes 27 


center 
rim 


0.84 ±0.08^ 
n ^0+0.04(3) 

u - oy -0.06 


0.16+ 



Cygnus Loop center 0.82t° m {4) 

northeast limb 0.22-0.30 (5) 



'mean density of the entire remnant. 

References. - (l)Rho & Borkowski (2002); (2)Harrus et al. (1997); (3)Enoguchi et al. 
(2002); (4)Miyata et al. (1998); (5)Miyata et al. (1994). 
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Table 8. Electron and ionization temperatures of mixed-morphology SNRs 



SNR name 


(region) 


kT c 
[keVl 

L J 


kT z 
[keVl 

L J 


Element b 


IC 443 a 

W49B a 

W28 a 

W44 

3C391 

Kes 27 


(center) 
(entire) 
(east) 
(center) 
(entire) 
(entire) 


n qq+ 04 

u - yy -0.05 
i 77 +0.04 
L ' '-0.03 

1 07+ - 43 

1 - u< -0.28 
n 7C+0-03 

0-53±8;S 
0.82t° : » 


1 4Q+0.09 

2 47+ - 21 

-0.20 

0.93lg;» 
0.68lg;« 
0.47lg;g 

o.83±g:i| 


S 

Ca 
Si 
Si 

Mg 
Si 



a Derived temperatures are for the high-temperature plasma. 

b This column shows elements used to estimate the ionization temperature of the plasma. 



Table 9. Electron and ionization temperatures of shell-like SNRs 



SNR name 


kT c 


log(ra c t) 


kT z 


Element reference 




[keV] 


[s cm' 3 ] 


[keV] 





Cassiopeia A 


o no+o.30 

°- uz -0.24 


10.831H 2 


1 n7+ - 02 

i - u ' -0.04 


s 


our work 


Kepler's SNR 


2.38 ±0.14 


10.84 ±0.03 


0.97 ±0.03 


s 


1 


Tycho's SNR a 


0.85 ±0.15 


11.0 ±0.1 


0.57 ±0.03 


Si 


2 


G352. 7-0.1 


2.ot :? 


ii.oi8:i 


l.UU_ 22 


s 


3 


Puppis A 


0.72 ±0.06 


10.47 ±0.07 


n oro+0.011 

U.OOO_Q QQg 


Mg 


4 


RCW 86 


0.430lffi 


10.71 ±0.05 


0.286 + °;^ 


Mg 


5 


Cygnus Loop b 


0.82i°;° 9 7 


10.35i°;° 6 7 


o.54±g:8i 


Si 


6 



References. — (l)Kinugasa & Tsunemi (1999); (2)Hwang & Gotthelf (1997); (3)Kinugasa 
et al. (1998); (4)Tamura (1994); (5)Rho et al. (2002); (6)Miyata et al. (1998). 

a For Tycho's SNR, both kT e and n c t are derived from Si and S line diagnostics, not from 
NEI model fittings. Details are described in Hwang & Gotthelf (1997). 

b isted kT e and n e t are those in the central region of the remnant, see Miyata et al. (1998). 



